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Abstract 

We investigate the phase structure of non-commutative scalar field theories and find 
evidence for ordered phases which break translation invariance. A self-consistent one- 
loop analysis indicates that the transition into these ordered phases is first order. 
The phase structure and the existence of scaling limits provides an alternative to the 
structure of counter-terms in determining the renormalizability of non- commutative 
field theories. On the basis of the existence of a critical point in the closely related 
planar theory, we argue that there are renormalizable interacting non-commutative 
scalar field theories in dimensions two and above. We exhibit this renormalization 
explicitly in the large N limit of a non-commutative O(N) vector model. 
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1 Introduction 



Non-commutative field theories have recently received a great deal of attention [1, 2, 
3, 4, 5, 6, 7, 8], stemming in part from the fact that they arise as low-energy descrip- 
tions of string backgrounds with anti-symmetric tensor fields [4]. Renormalizability of 
non-commutative theories remains an open question. 1 Standard approaches to demon- 
strating perturbative renormalizability (see for instance [5]) encounter difficulties be- 
cause of infrared singularities. These singularities are not associated with any massless 
propagating fields in the theory, but instead arise through loop effects. They can have 
dramatic physical consequences: for instance, in a theory whose classical action is that 
of a massive scalar with cubic interactions, the = vacuum becomes not just globally 
unstable (on account of an effective potential which is unbounded below), but locally 
unstable, as if the scalar had become tachyonic [6]. 

The main interest in this paper will be in scalar theories with 4 interactions. For 
the most part we will restrict attention to Euclidean signature and to even dimensions. 
Following [6], we will briefly review in Section 3.2 how non-planar one-loop graphs lead 
to a singularity in the one particle irreducible (1PI) two-point function: T^(p) — > oo as 
p — > 0. What this amounts to physically is long-range frustration: (0(:r)0(O)) oscillates 
in sign for large x. A natural expectation, given such a correlator, is that the usual 
Ising-type phase transition, to an ordered phase with (0) ^ 0, will be modified to a 
transition to a phase where (0) varies spatially. This is indeed what we will find in 
Section 4: more particularly, we will find a fluctuation-driven first order transition to 
a stripe phase, where only one momentum mode of the scalar field condenses. 2 In 
Section 4.6 we will consider more complicated ordered phases, where more than one 
momentum mode condenses. In the perturbative regime of Section 4, it turns out 
that stripe phases are favored; however, it appears that as couplings are increased, the 
system alternates between preferring the condensation of one or several momentum 
modes. 

The overall picture we will find for the phase diagram is a first order line terminating 
on one end at a Lifshitz point, where first order behavior merges back into the second 
order transition of the Ising model; and terminating on the other end at a critical 
point which arises in a planar version of the commutative theory. We review in Sec- 
tion 2 the relationship between phase structure and renormalizability, and argue via 

1 As this paper was nearing completion, we received [9], which argues that a four-dimensional non- 
commutative Wess-Zumino model is renormalizable. The argument hinges on controlling the infrared 
singularities which give rise to the interesting phase structure explored in this paper. Thus there is 
little overlap between [9] and our work. 

2 In [10] it was argued that a translationally invariant ordered phase with massless Goldstone bosons 
is impossible in continuum renormalized perturbation theory. This can be regarded as a hint of an 
exotic ordered phase such as the stripe phase that we find. 
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scaling that the existence of the critical point in the planar theory should imply the 
renormalizability of the non-commutative theory. 

The main method we use to establish the existence of phase transitions is a self- 
consistent Hartree treatment of one-loop graphs. This same method is our primary tool 
in demonstrating the validity of the scaling arguments that guarantee renormalizability. 
Thus, these arguments are airtight only for the large N limit of a non-commutative 
O(N) vector model, where the self-consistent one- loop Hartree treatment becomes 
exact. 3 More sophisticated methods are called for to decide the validity of scaling re- 
lations in more general non-commutative field theories. As explained in Section 5, it 
is possible to arrange the quartic couplings in the O(N) theory so that there are no 
divergences at all at leading order in N, independent of the dimension. Our renor- 
malizability arguments apply, however, for arbitrary quartic couplings in the O(N) 
theory. 

Other authors [13, 14, 15] have discussed finite temperature effects in non-commutative 
field theories, using the Matsubara formalism with periodic Euclidean time. The aim of 
this paper is rather different: we work with non-commutative field theories on uncom- 
pactified flat space, and varying "temperature" is regarded as equivalent to changing 
the bare mass. 

We conclude in Section 7 with a discussion of the relevance of our results to string 
theory and to quantum hall systems. 



We begin with an extremely brief recapitulation of the Wilsonian connection between 
the phase structure of a (cutoff) field theory and that of taking its continuum limit— 
which is the problem of renormalizability. As we are interested in this paper in the 
renormalization of scalar fields, we will recall the lore on commuting scalar fields. 

In the Wilsonian approach, we study the Euclidean field theory governed by the 
action 



as a statistical mechanics problem — i.e. we introduce a momentum cutoff A and mea- 
sure all dimensionful quantities in its units. This leaves us with a problem with one 
degree of freedom per dimensionless volume (A~ d in physical units) and dimensionless 

3 Strictly speaking, this is true only in the disordered phase and at critical points at the boundary 
of this phase. In the ordered phase the distinction between the self-consistent N = 1 problem and 
the large N limit could be important, see e.g. [11] and especially [12]. We expect to return to this 
question in the current context in the future. 



2 Phase Structure and Renormalizability 
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Figure 1: The phase diagram for the Ising scalar field theory, (1), in d > 1. 

couplings m 2 /A 2 and g 2 /A 4_d , commonly labelled r and u in the statistical mechan- 
ics/condensed matter literature (see for instance [16]). We next search for lines of 
continuous phase transitions, or critical surfaces, in the r, u plane; the critical surface 
in this problem (which is in the universality class of the Ising model) is sketched in Fig- 
ure 1. On this surface, the correlation length, in dimensionless units, £ (the "lattice" 
correlation length, literally so if the cutoff is implemented via discretization) diverges 
which is the sine qua non of taking a continuum limit. Having located a critical surface, 
three continuum limits are possible at each point on it — a massless limit obtained by 
sitting exactly on the critical surface, and two massive limits in which r(A) and u(A) 
are chosen to approach the critical surface from either phase as the cutoff is taken to 
infinity while keeping A/£(r,u) = M R fixed and equal to a renormalized mass. That 
such limits can be taken, requires scaling in the statistical mechanics. Also, the phe- 
nomenon of universality will imply that some domain of a critical surface will exhibit 
the same long distance correlations and hence give rise to the same continuum limit. 

In the above description we have used the language of phase structure. A more 
powerful account is that of renormalization group (RG) flows which we have sketched 
in d — 4 (Figure 2a) and d < 4 (Figure 2b). In the RG description we are interested 
in fixed points of infinite correlation length. In d = 4 we see that all critical theories 
flow into the gaussian fixed point (r, u) = (0, 0) whence the (strict) continuum limit is 
trivial, while in d < 4 the non-trivial Wilson-Fisher fixed point leads to an interacting 
continuum limit. This information is not available from a phase diagram alone. The 
fixed point analysis also implies scaling and hence guarantees a continuum limit. 

We belabor this point because in the balance of the paper, we will deal in phase 
diagrams and not RG flows. We will find critical points and will argue that these guar- 
antee the existence of continuum limits (renormalizability) of non-commuting scalar 
theories. In a specific large N theory, we will be able to show this explicitly, and we do 
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Figure 2: Renormalization group flows for the Ising scalar field theory in (a) d — 4 and 
(b) d < 4. 

not doubt that the claim is correct. Nevertheless, it does not have the full generality of 
an RG analysis and it would be nice to carry out such an analysis, even perturbatively 
as done for the commutative 4 theory in d = 4 by Polchinski [17]. 



3 One Loop Action for Non- Commutative Scalars 
3.1 Generalities 

Our starting point is the non-commutative scalar field theory specified by the action 



S 



I 



d d x 



1 2 ,2 , 9 2 



s*4 



T T ' 2 T 4 
where 0* 4 = * * * and the star product is defined as usual by 

(f*g)(x) = e^ 9 ^f(x)g(y) 



y=x 



(2) 



(3) 



The anti-symmetric matrix whose relation to the star product can be most simply 
expressed through 

[x^,x% = x» * x" - x" * x» = iQ^, (4) 



will for most of this paper be assumed to be of the form 

1 



0^ = 6 



-1 



<-d/2 



(5) 
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for even dimensions d. We will comment in Section 4.5 on the more complicated cases 
where Q^ u has unequal eigenvalues or d is odd. 
A useful extension of (3) is 



(/i * h * • • • * ft)(x) = e^«J e "^7i(*i)/ 2 (*2) • • • fe(x e ) 



(6) 



The result (6) is easiest to obtain via Fourier analysis, using the basic result 

where by definition p A q = Q^p^Pv (Condensed matter readers should note that this 
is the lowest Landau level algebra of density operators [18]). 

Like matrix multiplication, the star product is non-commutative. However, a product 
of exponentials e ipk ' x , can be reordered cyclically if the pk sum to zero. As a result it is 
necessary to specify the cyclic order of vertices in writing down Feynman rules, just as 
in large iV theories. It is well known [1] that planar amplitudes of the field theory (2) 
are independent of Miy (and hence are the same as in the commutative theory where 
= o), up to an overall phase. The effects of non-planar diagrams were first studied 
systematically in [6]. 

We will also consider two variants of (2), namely complex scalars, where 



/ 



d d x 



<90<90* +m 2 00* + ^0*0**0*0* + ^0*0*0**0* , (8) 



and the non-commutative O(N) vector model, 



1 1 q 2 q' 2 



S = J d d x 

The field theory (8) has been studied previously [19]. 



(9) 



3.2 One loop diagrams 

In this section we review the results of [6] which will be relevant for our calculations. 
The one loop corrections to the propagator of the scalar field split into planar and 
non-planar parts: 

planarW ^ J (2nYk 2 + m 2 
(2n) d k 2 + m 2 - 

Using the Schwinger parametrization 

I poo 

T^T-2= / dae~^ 2 \ (11) 
k l + m 2 Jo 
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one easily extracts 

rA d d k e ipAk 



{2n) d k 2 + m 



2 



00 da _„ m 2_P££ 1 



(An) d / 2 Jo ( 12 ) 



2-d 



= (27r) d/2 m 2 ^po]j+— j l^mypo p+ — 

where, following [6] , we have introduced an ultraviolet regulator A and defined a sym- 
metric product 

poq = -p lt & M/ G v x q x . (13) 

K u (x) denotes a modified Bessel function. 

The end result is a corrected propagator of the following form: 

T^(p) = p 2 + m 2 + 2<? 2 / 2 (0) + g 2 I 2 (p) . (14) 

The parallel results for the complex scalar and the O(N) vector model are 



complex scalar: T (2) (p) = p 2 + m 2 + (2g 2 + g' 2 )I 2 (0) + g' 2 I 2 (p) 
0(N) model: T (2) (p) = p 2 + m 2 + g 2 NI 2 (0) + g' 2 NI 2 (p) + 0(1) , 



(15) 



where for the 0(N) model we have only evaluated the graphs which contribute at 
leading order in large N. The salient property of r^ 2 ^(p) for our subsequent discussion 
is that it grows for small p. If we think of absorbing the planar one-loop contribution 
to r( 2 )(p) into the definition of mass (for instance, define a renormalized mass through 
M 2 = m 2 + 2g 2 I 2 (p 2 ) for the real scalar theory), and then removing the cutoff while 
holding M 2 fixed, then T^(p) is finite for finite p but goes to +00 as p — > 0. This 
means that the low-momentum modes are extremely stiff; thus they seem likely never 
to participate in a phase transition. If there is a phase transition, it should involve the 
momentum modes where T^ 2 \p) is smallest, since these are the ones most likely to be 
destabilized as M 2 becomes negative. The details of how this happens turn out to be 
somewhat intricate, and the next section is devoted to sorting them out. 



4 Phase Structure for Non-commutative scalars 
4.1 Generalities 

The non- commutative action (2) is characterized by three dimensionless parameters: 
m 2 /A 2 , g 2 /A 4 ~ d and 9 A 2 . Accordingly we need to establish a phase diagram in a three 



6 




Figure 3: The phase diagram for the non-commutative 4 theory at fixed g 2 , in d — 4. 
The solid lines show transitions for the case of a single real scalar. The point marked 
L is a Lifshitz point. The transition between uniform and non-uniform ordered phases 
is not accessible, in d — 4, via the techniques used in this paper. Instead, the dashed 
lines represent the transitions that arise in a self-consistent one-loop analysis. 



dimensional space. In what follows we will typically work at a fixed dimensionless 
coupling and take two dimensional cuts through parameter space instead. The resulting 
phase diagram is sketched in Figure 3, and the rest of this section is devoted to the 
arguments that give rise to it. 

4.2 The Planar Limit: 6A 2 = oo 

It has been noted previously that the large 6 limit picks out planar diagrams [1]. In 
the context of a cutoff field theory, this statement can be made precise. 

Consider the theory at m 2 / A 2 > so that all diagrams are infrared finite in addition 
to being ultraviolet finite on account of the cutoff. The perturbative expansion for the 
free energy F = — log Z is the sum of all connected vacuum diagrams (no external 
legs). Of these, the planar graphs involve no phase factors stemming from the non- 
commutativity while the non-planar graphs involve at least one. In the limit 8A 2 — > 
oo the latter vanish on account of the infinitely rapid oscillations of the integrand. 
Consequently, the high temperature expansion of F in the maximally non-commuting 
scalar theory is identical to that of the planar theory, defined as the sum over planar 
diagrams alone. 4 

4 For us, high temperature is synonymous with large positive m 2 , while low temperature means 
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This statement can be extended to selected classes of fields. The perturbative ex- 
pansion for the free energy in a field is, 



F[.J}= log Z[J] 



logZ[0]+W •••/ G c (p 1 , ....,p n )J(p 1 )....J(p n ) 

r, J PI JPn 



(16) 



n J Pi Jpn 



where G c (pi, ....,p n ), inclusive of the momentum conserving delta function, is the con- 
nected n-point Green function in momentum space. On account of the vanishing of 
the non-planar graphs alluded to earlier, 



for each n. For fields that are modulated only in one direction, the momenta p\ through 
p n are parallel whence the explicit phase factor for planar diagrams vanishes. Hence 
in the limit OA 2 — > oo, F[J\ is still given as a sum of planar diagrams alone. Note that 
for more complicated field configurations, the equivalence is no longer true. 

Thus far we have made statements in the high temperature phase. The planar the- 
ory, inheriting the Ising symmetry (<f) — > —<f)) of the full theory, will exhibit a twofold 
degenerate broken symmetry phase. This will be reached on traversing a continuous 
phase transition at a critical ml/ A 2 < which will be mean field in d > 4, mean 
field corrected by logarithms in d — 4, but likely in a different universality class from 
the standard Ising transition in d < 4 (see below). An immediate deduction from the 
equivalence of the infinite OA 2 and planar theories in a field is that we may analyt- 
ically continue their common free energy into the low temperature phase (bypassing 
the critical point), thereby establishing the equivalence of their low temperature ther- 
modynamics as well. Altogether, we may conclude that the infinite OA 2 theory has 
a critical point at the same m 2 /A 2 as the planar theory, which is exactly the same 
transition. This critical point will play a central role in achieving a continuum limit 
for the noncommutative theory. 

Finally we should note the well known result that the planar theory is also the 
N — > oo limit of a hermitian matrix model with an appropriately generalized version 
of (1) as its action. Via this route, there are exact results on the phase transition in 
d = 1 [20] 5 that show that the transition is different from that of the Ising model, 
and an approximate RG treatment that has yielded exponents in 2 < d < 4 as well 
[21, 22]. There does not appear to be a treatment of the broken symmetry phase in this 
formulation of the problem — at issue is how the Ising symmetry of the planar theory 

large negative m 2 . 

5 It is common to refer to the c = 1 matrix model as string theory in two dimensions, on account 
of the anomaly-induced dynamics for the Liouville field. From the point of view of the matrix model, 
however, d= 1. 



lim e ^<i PiApj G c (p 1 ,....,p n ;6A 2 ) 



Gf^( Pl ,....,p n ) 



(17) 
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is embedded in what appears prima facie to be a much larger symmetry group in the 
matrix model. 

4.3 First Order Transition at Large 6A 2 

Having established that there is a critical point at 8 A 2 = oo, we will now show that this 
terminates a line of fluctuation driven first order transitions at large 8 A 2 . That this 
should happen, is a fairly general expectation from the work of Brazovskii [23] when 
combined with the observation [6] that the one loop T^(p) has a minimum at non-zero 
p at large OA 2 (the latter condition is implicit in their analysis), and we will review 
this physics below. Following Brazovskii, we will employ the self-consistent Hartree 
approximation. This approximation is sensible in d > 4, but does not capture the 
planar theory transition in d = 2. Consequently, we will not treat that case in detail, 
although our later arguments on renormalizability will apply there as well. 

We will present two separate arguments for the first order transition. The first will 
involve an asymptotic construction of the solution to the self-consistent problem — this 
will follow the original analysis as closely as possible but with complications that we 
detail below. This will also require that we take a double limit in which OA 2 is taken 
large but g 2 is taken to zero — i.e. we will be working in the vicinity of the massless 
Gaussian theory. The second will be a more indirect argument which will be carried 
out at fixed g 2 and will consist of showing that the free energy of two solutions must 
cross in a certain region of parameter space. 

4.3.1 Take I 

We turn to the direct construction. To keep the discussion simple, we will continue to 
assume even Euclidean dimension d, and Q^ u with maximal rank and eigenvalues ±8. 
We will also omit inessential factors of order unity. 

A self-consistent treatment of the one-loop correction to the propagator leads to 

Having of maximal rank, with eigenvalues ±8, leads to a considerable simplifica- 
tion: r< 2 >(p) is a function only of p 2 . This can be shown inductively, order by order in 
g 2 . Suppose that T^(p) is SO(d — l) to £-th order in g 2 . The (£+l)-st order expression 
for r( 2 )(p) is then invariant on SO(d — 1) transformations of p^Q^, and so can only 
be a function of p op = 9 2 p 2 . Equation (18) can thus be re-expressed as 

r^H^+^w^^ * (19) 
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Figure 4: 1PI graphs with dressed propagators. 



Here m 2 is the bare mass, and 




d d k 1 



(20) 



is the mass corrected by the planar graph, Figure 4a. The non-commutativity does not 
affect the large p behavior, T^(p) ~ p 2 . Thus the integral in (20) must be regulated, 
for instance with an explicit cutoff as in Section 3.2. No other ultraviolet divergences 
will arise in the following computations, so we may in effect take the cutoff to infinity 
and treat M 2 as the finite quantity which we dial to produce a transition. This leads 
to an argument for renormalizability, as we shall explain in Section 6. 

The last term in (18) has a contribution from large momenta k which goes as l/p d ~ 2 
as p — > 0. This prevents condensation of very low-momentum modes of 0. Thus, if 
there is a phase transition at all, the ordered state must break translation invariance. 

From now on, let us restrict attention to d — 4. At the end of this section we will 
remark on the extension to other dimensions. Proceeding naively, we could solve (19) 
to the first non-trivial order in g by replacing T^ 2 \k) on the right hand side by k 2 + M 2 . 
Then as M 2 is decreased below — ci|, for some (easily computable) constant c\ of order 

unity, T^ 2 \p) becomes negative near p = p c ~ \fgjd, signalling a second order phase 
transition to an ordered state where some momentum modes of condense. Our aim in 
the rest of this subsection is to show that this naive result is altered in a self-consistent 
analysis to a rather more interesting conclusion: there is a first order transition to a 

7 /3 

stripe pattern at M 2 = — Ci| — Ci 3 -^, where C2 is another constant of order unity. 

First observe that (19) makes the second order transition impossible. If we were to 
suppose that T^ 2 \p c ) = for some M 2 and some p c , then the integral in (19) diverges. 6 
The possibility of a first order transition was first realized by Brazovskii [23], given a 
suitable effective Lagrangian, with a minimum in the inverse propagator at non-zero p c . 

6 To be more precise, the integral in (19) diverges if T^ 2 \p) <~ (p 2 — p 2 ) 2 for p near p c : that is, 
I^ 2 ) (p) cannot at the same time remain smooth and have a zero. To cement the conclusion that there 
are no second order phase transitions, one must check that it is also impossible for r^ 2 ^ (p) to have a 
cusp form such as \p 2 — p 2 \ 2v with v < 1. Such a form can eliminate the divergence in (19). But there 
is still a contradiction: dT^ 2 '(p)/dp diverges as p — > p Cl but an integral expression for dT^ 2 '(p)/dp 
remains finite. 
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In [23], such an effective Lagrangian was simply assumed, and a one- loop fluctuation 
analysis with this starting point was shown to lead to a first order transition. This 
same method was further developed in [24]. In our case, the basic Lagrangian does not 
have an appropriate inverse propagator with a minimum at nonzero p c — as we have 
seen, this only arises in the 1PI effective action after a one-loop analysis. One cannot 
feed T^ 2 \p) directly into a Brazovskiian analysis: all fluctuations are supposed to be 
incorporated into r( 2 ) (p) already. However, we shall see that a self-consistent treatment 
based on (19) reproduces the essential features of [23], and does support the conclusion 
that a first order transition takes place. 

Although r< 2 >(p) can never vanish, it should be possible to make its minimum as 
small as we please by decreasing M 2 sufficiently. We then have the approximate forms 

r^(p) ps p 2 for large p 

T {2 \p)^tl{p 2 -pi) 2 + r lovp^p c ^ 

Y^ip) ~ — for small p. 

[6p) 2 

Two regions of the integral in (19) will then make significant contributions: the large 
p region and the p ~ p c region. The full integral can be approximated as a sum of the 
contributions from these regions. The parameters p c and £ can then be determined 
self-consistently: for p ps p c , 

T^ 2 \p) ps p 2 + M 2 + g 2 / k 3 dk F >' -\- p n Fc> / 

\Jo k 2 9p 2 c Jo r + H(k 2 -p 2 c y 



,2 i n n iri'i 



p' + M z + J*— + 



(Op) 2 ^ 
Pc = ^ = M r = 2p 2 c + M 2 + 93/9 



(22) 

The last line of (22) follows from matching the previous line to the approximate form 

again suppressed inessential factors of order unity in 
(22), and we have assumed 3 < 1, so that <^ 1. This latter inequality will be 
essential to our further analysis. It turns out to be difficult to get the factors of order 
unity right, because £0 is only on the order of l/p c - the breadth of the minimum of 
T^(p), in the approximation scheme we have used, is comparable to its distance from 
the origin. 

So far we have done all calculations as expansions around the disordered phase, 
where (0) = 0. The putative ordered phase is a stripe pattern: (0) = 0o = Acosp c x. 
(Strictly speaking, we do not expect an exact cosp c x dependence, only a function of 
x with period 2n/p c . But because the momentum modes near p c make the dominant 
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(a) 



(b) 



(c) 



(d) 
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Figure 5: Additional 1PI graphs for the ordered phase. Dashed external lines indicate 
insertions of the background field O - 

contribution, neglecting higher harmonics should not change the qualitative picture. 
One should in principle consider sums of the form cosp c x + cos p c y, as well as more 
complicated superpositions; however, as we shall see in Section 4.6, the simple stripe 
solution is favored at small coupling). Writing = (fio + rj and expanding in fluctuations 
of r], one obtains 



We will treat perturbatively all terms with factors of both r\ and 4>q. This is justified 
if A is small, which is not completely obvious given that the transition is first order. 
But for small g, A is small very close to the phase transition — at least, this is an 
assumption which can easily be verified at the end of the computation. Classically, one 
would demand that the terms in (23) linear in t] should vanish: this corresponds to 
summing the graphs in Figure 5a and 5b. At the one-loop level, the graphs in Figure 5c 
and 5d also contribute, and the end result is 



where $ is the free energy per unit volume. Adding the graphs in Figures 5e and 5f 
to the self-consistent one-loop graphs that contributed to (18), one finds only a slight 
modification of the propagator in the disordered phase: 7 To (18) one must simply add 

7 Actually, the propagator is not diagonal in momentum space once the graphs in Figure 5e and 5f 
are added. Non-diagonal terms would have to be taken into account if we wanted to determine the 
exact form of the ordered phase; however they should not change the qualitative features of the phase 
transition. 




(23) 



(24) 
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the terms g 2 A 2 (l + |cosp A p c ). The first of these terms comes from the graph in 
Figure 5e, while the second is from the graph in Figure 5f. For momenta p on the order 
oi p c , the argument of the cosine is small (again assuming 9p 2 c is small), so for such 
momenta we may alter (22) to 

T^(p)^p 2 + M 2 + T f T2 + f^ + g 2 A 2 

{dp) 2 toVr (25) 

^ = T = \/f r = 2 Pc + M 2 + f+ T + g 2 A 2 . 

To summarize the situation so far: a self-consistent analysis of the graphs in Figures 4 
and 5 led to an expression for the minimum, r = r( 2 )(p c ), of the 1PI propagator T^(p); 
and an expression for the tadpole for the magnitude A of the condensate. These 
expressions are 

r = r + g ^ + ^ -- = A \r-'-A>), (26) 
where we have defined 

r = 2p 2 c + M 2 a = ^ Pc = \le- ( 2? ) 

The various factors of order unity which we have neglected can be absorbed by rescaling 
g, Pc and a. 8 None of these rescalings affects the conclusion that there is a first order 
transition. However, the rest of the analysis is somewhat delicate, and we must from 
now on be meticulous about factors. The calculation we are about to sketch appeared 
in [23] with some slight errors, which were corrected in [25], modulo a typo in the sign 
of the last term in their (2.13). 

A stable phase must have a vanishing tadpole. This happens either when A = (the 
disordered phase) or when r — \A 2 = (the ordered phase). The respective values of 
r, as well as the value of A in the ordered phase, can be determined from 

-r = 

(28) 



Note that the ordered phase only exists for r < — 3(a/2) 2 / 3 xs —1.89a 2 / 3 . In this range, 
there are two solutions for r Q (both positive): let us call these r ol and r o2 . The free 

8 In point of fact, it is difficult to compute some of these factors, since they arise from the integral 
equation (19). We thank D. Priour for consultations on the possibility of treating (19) numerically. 
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energy difference is conveniently calculated as 



A$ = $ - $ d = 




dA d$ 
dr dA 




(29) 



1 [ r 2 a +r 2 



This equation applies equally to the two solutions for r . By subtracting (29) with 
r o = r oi from the same equation with r Q = r o2 , one can verify that the larger solution 
for r leads to the lower free energy. One can also show, directly from (29), that the 
disordered phase is favored for —2.03a 2 / 3 < r < — 1.89a 2 / 3 , and the ordered phase is 
favored for r < -2.03a 2 / 3 . 9 

In terms of the original variables, and again in d = 4, the phase transition to an 
ordered phase occurs when M 2 = — Ci| — C2 2 ^p, for some constants C\ and c 2 of order 
unity. It is straightforward to show that c\ = l/ir and 9p 2 = g/2n. But to obtain an 
accurate value for c 2 requires numerics. The g 7 / 3 correction is a measure of the extent 
to which the system avoids the second order behavior expected from a naive one-loop 
analysis. Higher loop corrections are expected to shift the critical value of M 2 9 by 



4.3.2 Take II 

We now consider a second route to establishing a first order transition at large 6A 2 
which exploits the intimate connection between the phases of the noncommutative 
theory and their limits in the planar theory. The basic idea is this: As shown earlier, 
one cannot have a continuous transition once the disordered phase propagator has a 
minimum away from zero which is the case once 6A 2 is reduced from infinity. While 
this allows the high temperature solution to exist below the critical temperature of the 
planar theory, there must be a second solution that grows out of the ordered solution 
of the planar theory but is now modulated at a very long wavelength. If we take OA 2 
sufficiently large at fixed m 2 < m 2 , the ordered solution will win for it is the correct 
solution at 9 A 2 = oo. Hence there must be a first order line emanating from the planar 
theory critical point. 

To construct a proof on these lines it would appear that we could get away by com- 
paring free energies exactly at 9 A 2 = oo for continuity would give us a range of 9 A 2 

9 In principle one should be able to prove these assertions using the formula for solutions of a cubic. 
However, since the problem has only one parameter, r/a 2 / 3 , we have found it more convenient to 
proceed numerically. 



0(g 3 ). 
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over which the ordering would still hold. This almost works — while there is no problem 
with the ordered solution, which connects to the ordered solution of the planar theory, 
the delicate point is that the OA 2 — > oo limit of the disordered solution is not itself a 
solution of the planar theory self-consistency equation. Consequently we need to take 
the limit carefully, which can be done, modulo a weak and entirely plausible assump- 
tion, without actually solving the equations. (For the simpler Brazovskii problem, a 
similar strategy is entirely successful. This is detailed in the Appendix.) 

To begin, let us consider the limiting ordered free energy. This is the free energy of 
the planar theory, 

F ord = (-0 O + T O ) — 2 + - J j~y d log[* + g 0J (30) 

where O = (0) and is determined by m 2 via, 

d d k 



m2 + 4 ^ + 8g2 JWyWTM = - (31) 

In writing this compact form, we have used the relation £ = — m 2 + Ag 2 4> 2 valid for the 
fluctuation propagator in the ordered phase. 

We next turn to the disordered phase self-consistency equation, now written in terms 
of the self energy: 

, , n 2 f d d k 1 + \e lkhp 

s < ri ^7(2 ^ + ,J + s W - (32 > 

We will assume that the integrals are regulated in the ultraviolet, but the the precise 
choice of regulator will not be crucial. Consider the following assertions: 

(i) At any A and 6 A 2 , £(0) = §E(oo). This follows from the rapid oscillation of 
the phase factor as p — > oo. At increasingly large 8 A 2 , this fall in £(p) will happen 
increasingly rapidly. Hence if we examine F^(p) = p 2 + m 2 + S(p) it must develop a 
minimum away from p = 0. As already noted this minimum will prevent a continuous 
transition. 

(ii) For stability, T (2) (p) = k 2 + m 2 + E(p) > 0. 

(iii) It follows that £(oo) and S(0) are positive. 

(iv) The combination m 2 + must change sign between zero and infinity whenever 
m 2 < m 2 and 6A 2 is sufficiently large. The proof is by contradiction. At p = it is 
positive by (ii). Were it to remain bounded away from zero we could bound 

1 1 

< 17, (33) 



k 2 + m 2 + E(fc) k 2 



and deduce that E(oo) < — m 2 and thence that E(oo) + m 2 < m 2 — m 2 < which 
contradicts the hypothesis. 
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(v) The minimum value of m 2 + E(p) must vanish as 9 A 2 — > oo. This follows from (i) 
and (ii). The minimum value will be achieved (or arbitrarily closely approximated in 
the case of purely monotone decreasing E(p)) ever closer to p = in this limit. To 
avoid violating (ii), it must vanish. 

Thus far we have not made any assumptions. Now we need to make a mild assumption 
to proceed further. 

(vi) We assume that m 2 + E(oo) vanishes in the large 9 A 2 limit as well. If E(p) 
is monotone decreasing, this follows from (v). Otherwise is is strongly indicated by 
continuity from higher temperatures. For m 2 > m 2 the limit is a positive constant, 
being the renormalized mass of the planar theory, which vanishes at m 2 = m 2 . It 
seems highly unlikely that at any fixed 9 A 2 ^ oo the large momentum behavior will be 
non-monotone as a function of temperature whence the conclusion. In sum, we assume 
that E(p) goes pointwise to — m 2 . Note that this is not a solution of the limit of (32). 

(vii) The remaining task is to evaluate the free energy in this limit. The free energy of 
the disordered solution is given by, 

1 r d d k E(fc) 1 r d d k , ri2 , 

As the limit of the disordered solution is not itself a solution, we have to be a bit 
careful about evaluating F^ s as 9 A 2 — > 00. Setting E(p) + m 2 = in the logarithm is 
unproblematic. The remaining integral can be evaluated by breaking it up as 

r d d k m 2 + E(£;) r d d k m 2 

1 (2vr) d k 2 + m 2 + E(fc) ~ J (2vr) d k 2 + m 2 + £(fc) ( ^ 

wherein the first term vanishes at large 9 A 2 and the second can be evaluated via the 
self-consistency equation. Finally, 

(m 2 ) 2 1 f d d k , ri2l 

It is not hard to see that this is higher than the free energy of the ordered solution for 
m 2 < m 2 by expanding the former to leading order in 2 , which is what we set out to 
prove. 

4.4 Small 9 A 2 : Ising Transition and Lifshitz Point 

We now turn to the opposite limit, namely that in which 9 A 2 is sufficiently small — we 
will quantify that below. The analysis in this limit is straightforward. Exactly at 
9 A 2 = we recover the commuting theory which has the standard Ising critical point. 
At small, non-zero values of #A 2 we can expand the exponential phase factor in the 
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quartic vertex in powers of J2i<jPi ^Pj- The leading term is the quartic interaction 
of the commuting theory, and the additional terms are clearly irrelevant at its critical 
fixed point in d > 4. This statement is true below d = 4 as well in the e expansion, 
for the additional terms possess a momentum structure that is not generated by the 
interactions already present at the Wilson-Fisher fixed point. 10 

This RG statement can be understood more prosaically from the one-loop compu- 
tation reviewed in Section 3. On examining the minimum in T (2) (A;) we find that it 
continues to be at k = at small 6A 2 whence the phase transition would still be into 
the uniform Ising symmetry breaking state. As a bonus, we discover that there is a 
critical value (6A 2 ) C ~ 1/g at which the minimum starts to move away from k = and 
that initially it grows as the square root of the deviation in OA 2 . While these precise 
claims will be modified in an actual theory of this region (which the straight one loop 
computation is not, being as it is only valid for m 2 / A 2 > 0), the general feature that the 
Ising transition to a uniform phase will give way to a first order transition (expected 
by continuity with what happens at large 6A 2 ) to a modulated phase at a Lifshitz 
point should be correct. Unfortunately, our self consistent treatment is incapable of 
producing a Lifshitz point at a finite temperature (m 2 /A 2 > — oo) in d = 4 (which 
is the lower critical dimension in this approximation) and so the Ising transition and 
the first order line meet only at m 2 / A 2 = — oo. Further investigation of this region is 
clearly desirable. 

4.5 Odd dimensions 

The phase transition we have described is dimension-sensitive. In d = 2, there can be 
no long-range order: the stripe phase is unstable to infrared fluctuations. 11 In d = 3, 
with only 9 12 ^ 0, the story is slightly more complicated, since T^(p) is no longer a 
function only of p 2 . The momentum modes for which (p) is a minimum lie on a ring 
in the p±~P2 plane. One may nevertheless argue as before that a second order phase 
transition is impossible: the integral equation for T^^p) implies that it is smooth; but 
then if T^(p) is to have a zero along the ring, it must be a quadratic zero, and such 
a zero renders the integral infinite for all momenta. A weakly first order transition 
to a stripe phase is the expected behavior (weak because only logarithmic divergences 
make the transition first order). The main point that must be checked is that the stripe 
phase is not destroyed by infrared fluctuations. This we will do in section 4.6. 

Suppose we add one more commuting dimension: that is, we work in d = 4 but with 
only O 12 and O 21 nonzero. T^ 2 \p) should again have a minimum for p on a ring in the 

10 This assumes we continue the wedge structure in some fashion between dimensions four and two. 
11 This statement is essentially an application of the Coleman-Mermin- Wagner theorem; however we 
will check it explicitly in Section 4.6. 
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Pi~P2 plane. Since this ring is now codimension three, it is possible for r^(p) to go 
smoothly to zero in the self-consistent Hartree treatment. So a second order transition 
becomes possible. What happens with a generic is less obvious, but again a second 
order transition seems possible. One might be able to rule out a second order transition 
on the basis of an instability in the one-loop corrected quartic interactions. It would 
be interesting to work out this case in more detail. 

4.6 Ordering beyond Stripes 

Our analysis of the ordering in non-commutative scalar theories has thus far been 
restricted to stripe phases. Following Brazovskii, we would expect these to be the 
correct solution in the self-consistent Hartree approximation at weak coupling. However 
this ignores two orthogonal possibilities. The first is the option of going between the 
stripe phase and the disordered phase by analogs of the smectic and nematic phases in 
liquid crystal physics. 12 This is conceivable on symmetry grounds but as it is beyond 
the reach of our techniques, we are unable to say anything definite on this score other 
than to note that the destruction of long range stripe order in d = 2 (see below) would 
imply that any ordered phase would necessarily be more symmetric. 

The second possibility is that, at lower temperatures, the stripe phase might give 
way to one in which two (or more) different momentum modes of might condense, 
leading to a pattern of squares or rhombi. An interesting question here is whether the 
noncommutativity can influence the state selection in an interesting way. The purpose 
of this section is to argue that stripes are indeed preferred over rhombi for small 6p 2 c1 but 
that rhombi (or more complicated patterns) become favored as 6p 2 is increased. Since 
9p 2 is small when the coupling g 2 is small, the more interesting ordered phases can only 
arise at strong coupling. We will also remark on the first order phase transition to a 
triangular crystal that arises when a 3 interaction is present. The analysis will be very 
rough in that we will use a model lagrangian that will incorporate the soft modes and 
nonlinearities in the problem at the tree level instead of attempting a full fluctuation 
analysis. However, the conclusion that stripes give way to more complicated patterns 
as 6p 2 is increased should be robust, since it relies simply on the phases that emerge 
from the star product. 

As a preliminary, and to introduce our model, let us demonstrate that long range 
order for a stripe phase is impossible in two dimensions. This is a well-established 
result [27], so we will be brief. In the case of even dimension d, and assuming 6^ u 
has all eigenvalues equal to ±9, non-commutativity contributes to the analysis only 
by generating a minimum in the propagator for non-zero momentum p c . Thus we can 

12 We should note, however, that in the extensively studied problem of A-B diblock copolymer melts 
[26] there appears to be no evidence for intermediate phases. 
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work with an effective action of the Brazovskiian form [23]: 



^efT 



/A (i Kl 



(d 2 + pl)4>} 2 + + 



A variant for complex scalars is 

S eS = Jd d x ^i\(d 2 + p 2 M 2 + kM 2 + ■ 



(37) 



(38) 



It is assumed that k± > and k 4 > 0. When k 2 < 0, there is pattern formation. The 
complex case is simpler, since the classical minima of (38) are plane waves, <fi = Ae ip ' x , 
with \A\ = \J — K2/K4 and \p\ = p c . We will focus on the complex case at the outset, 
and return to real scalars, and to 3 interactions, near the end of this section. 

The precise claim about absence of long-range order in the model (38) is that at any 
finite temperature, in d < 3, long-range order is destroyed by infrared fluctuations. To 
see this, suppose = A(x)e iPcxl where A(x) is assumed to be slowly varying. Plugging 
this ansatz into (38), one finds 



Sefr = J d d x ^«ij(2ip c a B i + d 2 x ±)A\ 2 + k 2 \A\ 2 + ^k 4 \A\ 



(39) 



The massless modes are those where A = ^ — k 2 / K^e 10 ^. We will use x 1 - to indicate 
the directions perpendicular to the unit vector x 1 . To quadratic order in 9, (39) reduces 
to 

hplW) 2 + (dl9f 



S e S 



j d d x ( 


KlK 2 










K 4 






(40) 



Because the action is fourth-order in derivatives, Coleman-Mermin- Wagner arguments 
are stronger than usual: an estimate of the fluctuations gives 



(9 2 ) oc / 



d d q 



Aplql + q\ ' 



(41) 



which is infrared divergent for d < 3. 

The above argument is not substantially modified in d = 2 by non-commutativity 
because the infrared divergence in (41) refers to physics at far larger length scales than 
the scale of non-commutativity or of frustration. However, in d = 3, the premise of the 
model is wrong: non-commutativity can only exist for two of the three coordinates, 
say x 1 and x 2 . If a stripe phase forms in the x 2 direction, the dispersion relation (up 
to dimensionful constants) is u)(q) ~ qf + q\ + q\. The integral / d 3 q/u(q) is now finite 
in the infrared, so the stripe phase is indeed stable, and our remarks at the end of 
Section 4.3.1 were justified. 
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Let us now proceed to the comparison of energetics for stripes/rolls and rhombi, in 
a non-commutative version of the Brazovskii model, 13 

5cff = / d d x C e g 

t i i x (42) 

£ e ff = Ul | (<9 2 + p2)0| 2 + K2 |0| 2 + - KA( j) * * 0* * 0* + _ K4 ,0 * 0* * * 0* j , 

with /ti > 0, k 2 < 0, and K4 + K4/ > 0. For simplicity, we will assume that <fi is 
modulated at most in two directions, x 1 and x 2 , which are rotated onto one another by 
the action of 9^. More complicated situations can be imagined since we must assume 
d > 2 for the analysis to proceed at all; however, most of the interesting physics should 
refer to dimensions paired together by 9^ u . 

Clearly, there are still roll solutions, = Ae ip ' x , to the equations of motion. Rhombi 
would arise from an ansatz = A 1 e Wl ' x + A 2 e tp2 ' x with \pi\ = \p 2 \ = p c , but with p 1 
and P2 linearly independent. This ansatz is not a solution to the equation of motion, 
but if we can show that it has a lower energy than the roll solution, it is reasonable 
evidence that rhombic patterns form. Plugging the trial function into (42), we obtain 



(I0I 2 ) 


= \A, 


r+ 


\A 2 


2 


(0 * * 0* * 0*) 


= \A, 


| 4 + 


\A 2 


j 4 + 2(l + cospi Ap 2 )|^i^ 2 | 2 


(0 * 0* * * 0*) 


= \A, 


! 4 + 


\A 2 


j 4 + 4|A 1 A 2 | 2 



(Aff) = ^(l^xl 2 + \A 2 \ 2 ) + ^4 + MG^ll 2 + l^| 2 ) 2 

+ (/?4 cos pi Ap 2 + k 4 /)|Ai^4 2 | 2 , 



(43) 



where (U) indicates that U is to be averaged over space. 

Clearly, rolls are favored over rhombi when k# > | /C4 j . In order for rhombi to be 
favored over rolls without making k 4 + k% < (which would destabilize the theory 
altogether), we need k 4 > and p\ A p 2 in the vicinity of tt. Because 

Pi A p 2 = 9p\ sin angle(pi,p 2 ) , (44) 

we need 9p 2 ~ 9\ to be finite. 

Suppose we indeed arrange values of the «j where rhombi are favored over rolls. The 
simplest case is K4 > but — 0. The curious fact is that a square lattice appears 
almost never to be the preferred pattern: if 9p 2 c < tt/2, rolls are preferred; whereas if 
9p 2 > tt/2, then among the trial wave-functions studied so far the one with the lowest 
energy has angle(pi,p 2 ) = sin~ 1 (2^p 2 /7r). We have not found actual solutions to the 

13 We thank E. Witten for a conversation in which the simplest case of this argument was worked 
out for small Op^. 
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equations of motion with the point symmetries of a rhombic lattice. Because of the 4 
terms, any such solution would necessarily be a combination of infinitely many plane 
waves with wave-vectors on the reciprocal lattice. It is straightforward though tedious 
work to optimize a variational ansatz which is the sum of finitely many plane waves. 
However there is not much point in going through the exercise because the form of (38) 
is only intended to capture the behavior of momenta close to \p\ = p c . 

Finally, let us turn to the case of a real scalar. We will continue to use the approach 
of a trial wave-function composed of plane waves with all momenta on the ring \p\ = p c . 
Plugging the ansatz 

= A ie ipi - X + A 2 e ip2 - X + c.c. , (45) 

into the effective action 



S cS = J d d xC cS 



1 



2 



1 2 1 



(46) 



2 ^ 4 



with Ki > 0, k 2 < 0, and k 4 > 0, one obtains the following: 



l ) = 2(\A 1 \ 2 + \A 



2 



(0*0*0*0) = Q(\A 1 \ k + |A 2 | 4 ) + (16 + 8 cos pi Ap 2 )\A 1 A 2 \ 2 

(C eS ) = k 2 (\A ± \ 2 + \A 2 \ 2 ) + -« 4 (|A| 2 + \M 2 ) 2 + «4(1 + 2cos Pl A p 2 )\A 1 A 2 \ 2 . 

(47) 

The sign on the last term can be negative, and the conclusion is that on the space 
of trial wave-functions that we have examined, rhombi are preferred over stripes if 
Pi Ap 2 > 27r/3 is finite. This last inequality is possible when dp\ > 2n/3. In a similar 
way, one can consider the condensation of three momentum modes and show that for 
9pl in a neighborhood of 2ir/ y^3, a hexagonal pattern is preferred over both rhombi and 
stripes. Again, although these trial wave-functions are not solutions to the equations 
of motion, it is reasonable to expect solutions with the same symmetries to exist and to 
compete with one another in a similar way. We have considered condensing momentum 
modes which all lie in a single plane on which 9^ v = 9e^ v . In d = 4 and higher, there are 
more complicated possibilities where momentum modes condense in many directions. 
We will not attempt to classify all the possible ordered phases. 
One can also consider adding a cubic term to (46): 

S e s = / d d x C c g 

1 r 12 l l 1 ^ 

£ c fr = [(d 2 + pD^I + 2 K2< ^ + 3 Ks ^ ~ k( f ) ~ k( f )+ * * * . 

As explained in [6], a 0* 3 interaction contributes a term —I 2 (p) to T^(p); thus if the 
cubic interaction is strong as compared to the quartic interaction, T^ 2 ^ (p) is monotonic, 
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and there is no reason to think that there are spatially non-uniform phases at all. For 
the purposes of this discussion, let us assume then that the cubic interaction is weak 
but nonzero. 

Given that the leading nonlinearity is cubic, the natural expectation in two di- 
mensions would be that commutative versions of (48) exhibit a first order transition to 
hexagonal crystals as n 2 is lowered. This persists in the presence of non-commutativity, 
as the following computation shows. The ansatz for is 

= A ie ipi - X + A 2 e ip2 ' x + A 3 e- i(pi+P2) - x + c.c. , (49) 

with \pi\ = \p 2 \ = \pi + P2I = Pc- In two dimensions, the only way that this can be 
arranged is for p 1 , p 2 , and —p\ — p 2 to point to the vertices of an equilateral triangle 
centered on the origin. The relevant spatial averages are 



2 > = 2£|A, 



|2 



(0 * * 0) = 3 cos ^-^{A 1 A 2 A 3 + A\A* 2 Al) 
(0*0*0*0) = 6 (£|^| 2 ) +4(l + 2cospi Ap 2 )Y,\ A i A 



(50) 



(£ cff > = K 2 J2\Ai\ 2 + k 3 cos ———^——( y AiA 2 A 3 + A\A* 2 Al) 

+ (El^l 2 ) +«4(1 + 2cos Pl Ap^YlAAjl 2 . 

As long as the cubic term is present, a first order transition to a hexagonal lattice is 
the expected behavior. To be more precise, the expectation is that with k 2 sufficiently 
small but positive, there is a minimum where (£ e s) is negative, and where all three A4 
are nonzero. To see that this must be so, set A { = An,i with X^n 2 = 1. Then 

(£ eff ) = k 2 A 2 + k 3 A 3 + k 4 A* 

K 2 = K 2 

k 3 = 2k 3 cos — - — n x n 2 n 3 (51) 



K 4 = K4 



3 

- + (1 + 2 cospi A p 2 ) n l n ) 



One may easily show that £4 > 0, so there is no runaway behavior. It is straightforward 
to show that (C c g) attains negative values precisely if ^^2^4 < 2|k 3 |. This is the 
desired result: even assuming a small cubic interaction, the first order transition at 
finite positive k 2 swamps the would-be second order behavior at k 2 = 0. The only 
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exception is when cos Ei/^i — q : a t this point, the cubic term vanishes. A solution 
with | Ai| = \A 2 1 = | A3 1 is still preferred, but the phase of A1A2A3 is a flat direction. 
On general Landau-Ginzburg theory grounds one might expect a second order point at 
cos 2i^£2 _ q se p ara ting the two first order lines corresponding, respectively, to phase 
locking of the three plane waves with AiA 2 A 3 positive or negative; however it seems 
likely that fluctuations again drive the transition first order. In higher dimensions the 
story is again more complicated. The usual expectation is a lattice with the maximal 
number of equilateral triangles; in the presence of non-commutativity these lattices are 
likely to have deformations and preferred orientations. We leave an investigation of 
the possibilities for future work. 



5 The non-commutative 0(N) vector model 

Naively, there are two reasons to think that quantum field theories defined on non- 
commutative spaces are under better control perturbatively than their commutative 
counterparts. First, non-planar loop diagrams include an oscillatory factor in the 
integrand which generically cures ultraviolet divergences [6]. Second, compactifying 
the non-commutative position space also entails a compactification in momentum space 
[28]: this is one of several manifestations of the interplay between ultraviolet and 
infrared effects. But there are also reasons to think that such quantum field theories 
are not under good perturbative control. First and perhaps most seriously, it has not 
been shown that the special form of the tree level lagrangian (polynomial in derivatives 
and fields, but with all ordinary products replaced by star products) is preserved by 
quantum corrections beyond one loop. Second, planar divergences are just as bad as 
for commutative field theories. 

The latter two difficulties can be resolved for special quantum field theories. Consider 
the O(N) vector model: 



-/ 



d d x 



1 1 q 2 q' 2 



, (52) 



which, in the large N limit, is dominated by bubble graphs in its disordered phase. If 
g 2 = 0, the loop integrals in these graphs converge even without an ultraviolet cutoff, 
due to the oscillatory factor in the loop integrands. Thus there are no counterterms 
in the large N limit: the theory is perfectly finite. At subleading orders in N, coun- 
terterms do arise, and all the usual questions arise regarding whether the special form 
of the lagrangian is preserved in the quantum theory. However it can perhaps be 
regarded as progress toward deciding the issue of perturbative renormalizability that 
these difficulties can be suppressed by 1/N in an appropriate 't Hooftian limit. 
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Returning to the bubble sum for (52) with g 2 = 0, this can be reduced to an integral 
equation — the 1PI two-point function is given implicitly by 

^^-'^ilw' (53) 

In other words, the self-consistent Hartree approximation in the disordered phase is 
exact at large N. The existence of an ordered phase, following our earlier caveat, 
remains an open question. 

It may be that (52) with g 2 = is a special case of a more general strategy for 
generating quantum field theories with improved convergence properties in a special 
limit — the main ingredient being two conflicting notions of planarity, one from the 
index structure of the fields, and one from the star product. However, we are not 
guaranteed to obtain a finite quantum field theory (even in special limits) through 
this trick, as there may be graphs which are planar in both senses. For example, the 
lagrangian 

C = \ + \nf tr e + £m*WM^ * •** * * •„.. , (54) 

where $ is an JV x iV hermitian matrix, specifies a quantum field theory with two 
conflicting notions of planarity; but there are still certain planar diagrams diverge, 
for instance the one-loop correction to the quartic interaction vertex. This divergence 
requires a counterterm of the form tr($ * $ * $ * $), so we learn that the special form 
of the lagrangian (54) is not preserved by quantum corrections. 

6 Continuum limits 

We turn now the question of taking a continuum limit, i.e. the question of renor- 
malizability. We will examine this in the disordered phase of the theory within our 
self-consistent Hartree approach. Formally, we will consider the renormalization of the 
infinite N theory, (52), with g 2 N and g' 2 N finite. In this approximation, the only 
non-trivial correlation function that enters the theory is Y^ 2 \p) but on account of the 
noncommutativity it has significant momentum dependence that makes the procedure 
somewhat more complicated than it might seem at first sight. Nevertheless, the claim 
is that there are non-trivial continuum limits, and therefore interacting renormalizable 
field theories, in any even dimension. It was conjectured in [6] that non-commutative 
4 theory should be renormalizable in d = 4, despite the infrared divergences. The 
results of this section do not amount to a demonstration of this claim, because we 
persist in working at infinite N; however we hope that an extension to finite N may be 
possible. We should note though, that our Wilsonian attempt to renormalize by means 
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of the planar theory critical point is closely connected to the the planar subtraction 
algorithm suggested in [6]. 

As a warmup recall the problem of the commuting, infinite N, O(N) vector model. 
Here, 

Y {2 \p) =p 2 + m 2 + T, (55) 

where 

2 r \ Jdh. i 

' ? (56) 



{2ix) d k 2 + m 2 + E ' 

where we have indicated the regulation of the integral explicitly. In 2 < d < 4, this 
theory has a critical point m 2 = m 2 (A) at which the renormalized mass m 2 + E vanishes. 
Explicitly, 

2 (? 2 / A rf d A; 1 

mc = -yi (27pF' (57) 

but that is not central. 

From the existence of the critical point it follows that we may choose m 2 (A) in order 
that m 2 (A) + X(g 2 , A) = M 2 is held fixed, 14 whereupon we may take A to infinity and 
obtain a renormalized propagator T^\p) = p 2 + M 2 . In this example the end product 
of this mass renormalization is trivial, but the underlying lattice problem is not. A 
non-trivial correlation length exponent v is hidden in the relation between the bare 
mass and the renormalized mass. 

Turning now to the noncommuting infinite N theory, we can divide S(p) into the 
parts coming from the planar and non-planar diagrams and rewrite the self consistency 
equation in two parts: 

M 2 = m 2 + 2g 2 N T f\ 10 * - (58) 
y J {2Tx) d k 2 + M 2 + H nv {k) y 1 

and 

Note that M 2 — m 2 = E(oo) introduced earlier. 

Now, the planar theory has a critical point at m 2 = m 2 and #A 2 = oo where M 2 
vanishes. Keeping a fixed dimensionful 9 then sends ^A 2 to infinity automatically, while 
the critical point enables us to choose m 2 (A; 9) such that M 2 is held fixed as A — > oo. 

In this limit we are able to remove the cutoff from (59) as well which still defines 
a finite E np (p) and thereby a finite renormalized two point function T^\p) = p 2 + 

This is the point in the argument where the existence of a critical point in the bare theory is 
crucial. As we know that we can tunc the bare mass to get M 2 to vanish at any value of the cutoff, 
it follows that we can hold it fixed at a specified value as well. 
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M 2 + £ np (p). Note that by dimensional analysis E np has the scaling form £ np = 
M 2 f(p/M,9M 2 ). Note also, that £ np ~ (gp) V 2 is divergent in the infrared. 15 

In this fashion we see that the existence of the planar theory critical point does 
allow a continuum limit to be taken at iV = oo. There is one more limiting theory in 
this case, the massless purely planar theory obtained by sitting exactly at the critical 
point. (If the N = oo theory has a broken symmetry phase then a third continuum 
limit would be feasible from within that phase.) We expect that the planar limit does 
not change much between N = oo and N = 1 in d > 4. That suggests that more 
progress could be made in establishing renormalizability via an 1/N expansion — the 
same logic underlies the comparable demonstration for commutative scalar theories. 
However, we have not examined this question in any detail. 

Note also the interesting feature that the continuum limit is nontrivial in the infrared 
while the ultraviolet behavior is free. The latter is consistent with the triviality of 
commuting scalar theories in d > 4 although in this case it cannot be distinguished 
from the vanishing of the anomalous dimension of the scalar field at iV = oo. At any 
rate, the point is that the continuum limit will be nontrivial even above d = 4 (say in 
d — 6) for finite noncommutativity is, in a sense, a relevant perturbation at the planar 
theory fixed point. This is in contrast to the situation in the commuting case where 
only the free massive theory is possible. 

Finally, we should note that a different set of continuum limits is possible near the 
Lifshitz point in the N — 1 theory. These would entail keeping OA 2 finite as A — > oo 
and hence a vanishing 9. We have not studied these, but our large N approach could 
be extended above d — 4 should interest in these theories be warranted. 

7 Concluding Remarks 

The purpose of this paper has been to investigate unusual phase structure in the 
simplest interacting non-commutative field theory, namely 4 theory. We expect the 
existence of stripe phases to be quite common in non-commutative theories, the reason 
being that r^ 2 ^(p) typically has singular behavior as p — > when the cutoff is removed. 
If T^(p) — > — oo as p — > 0, the theory is sick in the sense that for no value of bare 
masses have we found a stable vacuum. If T^(p) — > +oo, then some version of the 
arguments of Section 4 should establish a first order transition to a stripe phase. 

It is natural to ask, what manifestation might this transition find in string theory? 
At present, we have no definite answer; but let us remark on one obvious venue where 
such a transition might be expected to arise. D-branes in bosonic string theory can, 
in a rough approximation, be thought of as classical lumps of an open string tachyon 

15 Note that this is the same integral equation that arose in (53), albeit with M 2 replacing m 2 . 
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field T (see for example [29]). The potential for the tachyon field is cubic, so non- 
commutativity (in the form of a -B^-field) merely destabilizes the T = vacuum 
further. Unstable D-branes in type II superstrings, however, have a quartic potential 
for the tachyon field, and it is possible that a stripe phase will arise for sufficiently large 
B^y. If a stripe phase indeed exists, it would be quite a peculiar string background, 
not readily comprehensible in classical terms. Classically speaking, the stripe phase 
would amount to an alternating arrangement of stable D-branes and their anti-branes, 
parallel to one another. Such an arrangement seems obviously unstable toward the 
branes and anti-branes collapsing into one another. There are however some caveats: 
first, in working with the non-commutative field theory as an approximation to the 
string dynamics, we are by assumption working in a limit where the closed string 
interactions are suppressed. So the collapse of branes into anti-branes could have a 
much longer time scale than the transition to a stripe phase from a T = phase. 
Second, working directly with a field theory is suspect because the cutoff (l/y/a') is 
on the same order as the tachyon mass. A full string theory computation, with finite 
a', seems to be the only wholly reliable approach. 

We should note that our description of stripe phases may not incorporate an impor- 
tant piece of the physics of unstable D-branes: namely that the open string degrees 
of freedom are believed to be confined when the tachyon field has condensed. The 
mechanism for this confinement is not yet wholly clear (see [30, 31] for two interesting 
proposals), but it could play a role in the stability of various phases. 

A relatively well-understood aspect of unstable D-branes in a strong background B^ v 
field is condensation into a semiclassical configuration known as the non-commutative 
soliton [7]. It has been argued [8] that this configuration represents the collapse of a 
bosonic Dp-brane into a D(p-2)-brane; see also [32] for closely related work. The role 
of these objects in our quantum considerations is not entirely clear. For one thing, the 
stable solitons/instantons at infinite 9 found in [7] by leaving out the gradient term in 
the action, have exactly zero action in our case which makes it imperative to keep the 
gradient term in the analysis. Should such an improved computation be feasible, we 
would still suspect that in the large N limit these objects will be unstable to unwinding 
via other directions in field space and hence do not play a role. 

In the N = 1 case, the solitons/instantons could exist. As instantons we would 
have to consider whether they invalidate our conclusions on the nature of the ordered 
phases. This seems unlikely, at least at weak coupling where the wavelength of our 
stripe phases is much larger, by a factor of 1/ y/g, than the size of the instantons in the 
classical analysis. More generally, they are local distortions of the condensate, but do 
not (unlike, say vortices in an XY model) affect the ordering at large distances. Hence 
we expect that they will renormalize the properties of the ordered phases but not lead 
to any singular effects. 
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In contrast, actual solitons in a quantum theory in odd dimensions would rely on the 
existence of an ordered phase. In the classical analysis, this is assumed to be a uniform 
condensate but as we have seen this is no longer the case in the quantum theory. 
Again, at weak coupling, it is plausible that the solitons are essentially undistorted on 
account of their size, but their dynamics would clearly be affected by motion in an 
inhomogeneous background. 

It is perhaps worth recalling some features of quantum Hall physics that may find a 
formulation as noncommutative field theories. The first is the physics of the v — 1/2 
state which has an elegant microscopic interpretation in terms of dipoles 16 with evident 
parallels to the discussion of Bigatti and Susskind [3]. We note that stripe phases do 
occur in quantum Hall systems at half filling (in high Landau levels) [34, 35] but 
hasten to add that any connection between them, the dipolar v = 1/2 theory and our 
noncommutative results is purely speculative at this point. 

An equally speculative connection is to a widely studied model of the integer quan- 
tum Hall effect, namely Pruisken's non-linear sigma model. It has several incarnations, 
but the simplest has a lagrangian of the form 

£ = ^2 tr(9T) 2 + i<) tr (T[d x T, d y T\) , (60) 

where T takes values in the coset U(2N) / (U(N) x U (N)). The second term is topologi- 
cal. The model arises from a replica treatment of disorder in a theory of non-interacting 
electrons in a magnetic field. In the end, physical quantities must be computed in a 
iV — > limit. In the original formulation, the RG flow in the (a xy , cr xx ) plane was ar- 
gued to have fixed points at a xx = and a xy G Z and at a xx = 1/2 and a xy G Z + 1/2. 17 
The first set of fixed points represent the quantum Hall plateaux, and the second set 
control the critical behavior of transitions between plateaux. The behavior a xx — > in 
the infrared is the analog of confinement in this model, and the special fixed points at 
non-zero a xx are believed to arise from instanton effects, similar to 't Hooft's treatment 
of deconfinement in QCD at 9 = n. 

It is suggestive that precisely the U(2N)/(U(N) x U(N)) coset arises as the vacuum 
manifold in the condensation of unstable D-branes in type II string theory. However 
it is difficult to see how the topological term tr (T[d x T, d y T\) can arise in the string 
effective action. This term is the leading order term in a derivative expansion of 
tr (T* 3 — T 3 ); however, from a type II string theory point of view, neither the mixing 
of star products with ordinary products inside a single trace, nor the cubic form of 
this potential, is expected. Nevertheless, strictly from a field theory point of view, it 

16 For a recent discussion including an invocation of non-commutative geometry, sec [33]. 
17 However, Zirnbauer [36] has argued persuasively that the fixed point theories have instead a 
Wess-Zumino-Witten form. 
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would be interesting to ask whether replacing the topological term with tr (T* 3 — T 3 ) 
preserves the universality class. We hope to return to these issues in the future. 
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Appendix: The Brazovskii Transition by Non-Bra- 
zovskiian Means 

Here we will give a version of the argument used in Section 4.3.2 for the transition 
studied originally by Brazovskii [23]. The relevant action is (37), which we rewrite as 

S cS = J d d x C cS 

The traditional route to showing that there is a first order transition is a weak coupling 
analysis. Here we will work at fixed coupling but instead take p c to be the "small 
parameter" . For this to work we will need that the theory at p c = possess a continuous 
transition to a broken symmetry phase, else nothing is gained. In the self-consistent 
Hartree approximation of Brazovskii for the high temperature phase, 

g 2 r d d k 1 
~ ~2 J (2vr) d n(k 2 - pi) 2 + m 2 + E ^ ' 

this requires that we take d > 4. 18 In these dimensions, the critical point of the uniform 
(p c = 0) theory is at w 2 = -£ J 

Now consider p c > at fixed g 2 . By the standard argument sketched in Section 4.3.1, 
there is always a disordered solution at any m 2 /A 2 > — oo, with a free energy (written 
in terms of self-consistent parameters) 

F = + s / logW * 2 " p ' )2 + m2 + s| • (63) 



18 Alternatively, we can modify the Brazovskii action by replacing (fc 2 — p 2 ) 2 by (k — p c ) 2 which 
preserves the feature of a codimension one surface of soft modes while retaining the connection to a 
critical point in d > 2. 
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Figure 6: A first order line ending at a second order point. 



whence a continuous transition out of the high temperature phase is impossible. 

Next, we follow the disordered solution as we take p c to zero. For m 2 > ml the 
solution smoothly connects to the disordered solution of the uniform theory with £ + 
m 2 > 0. But for m 2 < m 2 c this is not possible as there is only an ordered solution. So 
in the vicinity of the p c = line we must have two solutions, where the second evolves 
out of the broken symmetry solution of the uniform theory as p c is tuned away from 
zero. 

The remaining task is to show that for m 2 < m 2 c there is always a value of p c 
below which the ordered solution wins — this then shows, as in the example of the 
noncommutative theory in Section 4, that there is a first order line in the m 2 ,p c plane 
(see Figure 6). For this again it suffices to make the comparison in the limit p c — > 
on the grounds that free energies are continuous. 

The disordered solution is readily shown to require that m 2 + £ -> as p c — > 
whenever m 2 < m 2 c . Hence its free energy has the limit, 




(64) 



The free energy of the ordered solution at p c = is 




log^A; 4 + y -<Pl] 



(65) 



where the parameters satisfy, 




(66) 



and 




(67) 
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An explicit comparison, easily done perturbatively in (p 2 01 shows that the ordered 
solution has lower energy, which is what we set out to show. 
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